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ABSTRACT 

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based 
on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 
to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency 
coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual 
shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were 
fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral 
indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an 
emphasis on proper, physical modelling of the synchrotron bump using multiple components. Planck ERCSC data also suggest that the original 
accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 
2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data 
indicate that gamma-ray emission must originate in the same shocks that produce the radio emission. 

Key words, galaxies: active - BL Lacertae objects: general - quasars: general - radiation mechanisms: non-thermal - radio continuum: galaxies 


1. Introduction 

This paper is part of the first series of publications based 
on measurements made with the Planck} satellite. Planck 


1 Planck (http://www.esa.int/planck) is a project of the European 

* Corresponding author: A. Lahteenmaki, alien@kurp.hut.fi Space Agency (BSA) with instruments provided by two scientific con- 


1 




Planck Collaboration: SEDs and radio spectra of northern AGN 


(Tauber etal. 2010; Planck Collaboration 2011a) is the third- 
generation space mission to measure the anisotropy of the cos- 
mic microwave background (CMB). It observes the sky in nine 
frequency bands covering 30-857 GHz with high sensitivity 
and angular resolution from 31' to 5'. The Low Frequency 
Instrument (LFI; Mandolesi et al. 2010; Bersanelli et al. 2010; 
Mennella et al. 201 1) covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HFI; 
Lamarre et al. 2010; Planck HFI Core Team 201 la) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled to 0.1 K. Polarization is measured in all but the 
highest two bands (Leahy etal. 2010; Rossetetal. 2010). A 
combination of radiative cooling and three mechanical cool- 
ers produces the temperatures needed for the detectors and op- 
tics (Planck Collaboration 2011b). Two data processing centres 
(DPCs) check and calibrate the data and make maps of the sky 
(Planck HFI Core Team 2011b; Zaccheietal. 2011). Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for Galactic and extragalactic astrophysics 
as well as cosmology. 

The paper uses data from the Planck Early Release Compact 
Source Catalogue (ERCSC; Planck Collaboration 2011c). The 
ERCSC provides positions and flux densities of compact sources 
found in each of the nine Planck frequency maps. The flux densi- 
ties are calculated using aperture photometry, with careful mod- 
elling of Planck's elliptical beams. The ERCSC includes data 
from the first all-sky survey, taken between 13 August 2009 and 
6 June 2010. This unique dataset offers the first glimpse of the 
previously unmapped millimetre and sub-millimetre sky. It is 
used here to create spectral energy distributions of 104 radio- 
bright, northern active galactic nuclei (AGN), with the most 
complete coverage in the radio to sub-millimetre frequencies to 
date. 

Radio-loud AGN host jets of relativistic matter emanat- 
ing symmetrically from the core. These jets produce copious 
amounts of non-thermal radiation, which dominates the spec- 
tral energy distributions (SEDs) of such sources compared with 
any thermal emission from the nucleus, i.e., the accretion disk. 
The SEDs typically consist of two broad-band bumps, the one 
at lower frequencies attributed to synchrotron radiation, and the 
other at higher frequencies attributed to inverse Compton (IC) 
radiation. The peak frequencies of the two bumps vary from 
one object to another. The peak of the synchrotron component 
can be between the infrared and X-ray domains, and the IC 
peak can range from MeV to GeV energies. The sequence in 
the peak frequencies of the emitted energy, and the factors that 
create it, have been a hot topic in blazar research for more 
than a decade (e.g., Fossatietal. 1997; Ghisellini et al 1998; 
Padovani 2007; Ghisellini & Tavecchio 2008; Nieppola et al. 
2008; Ghisellini & Tavecchio 2008; Sambruna et al. 2010). 

Research on blazar SEDs has concentrated on two ap- 
proaches. The first is fitting the SED with a simple function, 
usually a second or third order polynomial, to obtain the pivotal 
peak frequencies and luminosities of the radiation components in 
a straightforward manner. This approach is typically used when 
studying large samples (Fossati et al. 1998; Nieppola et al. 2006; 
Sambruna et al. 2006; Nieppola et al, 2008; Abdo et al. 2010a). 
The second approach is detailed modelling of the SED, starting 
with the definition of initial parameters such as electron energy, 
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magnetic field intensity, and Doppler factor. This method is more 
time-consuming and is used mostly for individual sources (e.g., 
Acciari et al. 2010; Collmar et al. 2010). The standard model is 
a leptonic, homogeneous, one-zone model, where the emission 
originates in a single component (for a review of the blazar emis- 
sion models, see Bottcher 2010). One-zone models are useful as 
first-order approximations, but, in reality, AGN jets are rarely, if 
ever, dominated by a single source of radiation. The material in 
the jets flows through shocks in the jet, which locally enhance 
the radiation (the “shock-in-jet” model; Marscher & Gear 1985; 
Valtaoja et al. 1992). There can be several of these shocks in the 
jet simultaneously, and adding these to the emission of the quies- 
cent jet, we have several radiation components. Therefore, ide- 
ally, the SEDs and radio spectra should be modelled with more 
than one component. Such modelling is also necessary for the 
proper identification of the high frequency (IC) emission sites. 

In this paper we use the nearly complete SEDs and well- 
covered radio spectra provided by Planck to look for signs of 
these multiple components contributing to the total radiation that 
we observe. The paper is structured as follows. In Sects. 2 and 3 
we introduce our sample and summarize the multifrcquency data 
used in our study. The spectral energy distributions and radio 
spectra are discussed in Sects. 4 and 5, and their modelling is 
described in Sect. 6. In Sect. 7 we discuss the implications of our 
results for understanding the acceleration mechanisms in blazar 
jets, and in Sect. 8 we summarize our conclusions. Throughout 
the paper we adopt the sign convention for spectral index, a: 
S y a v“. The errors of numerical values marked with a plus- 
minus sign correspond to one standard deviation. 

2. The sample 

The complete sample presented in this paper consists of 104 
northern and equatorial radio-loud AGN. It includes all AGN 
with declination > - 10° that have a measured average radio flux 
density at 37 GHz exceeding 1 Jy. Most of the sample sources 
have been monitored at Metsahovi Radio Observatory for many 
years, and the brightest sources have been observed for up to 
30 years. The sample can be divided into subclasses as fol- 
lows: 33 high-polarization quasars (HPQs), 21 low-polarization 
quasars (LPQs), 21 BL Lacertae objects (BLOs), 19 normal 
quasars (QSOs), 9 radio galaxies (GALs), and one unclassi- 
fied source (J1 849 15+67064). (See, for example, Hovatta et al. 
(2008, 2009) for additional information on the classification.) By 
high-polarized quasars wc mean objects which have a measured 
optical polarization > 3% at some epoch, while low-polarized 
quasars have a polarization <3%. Normal quasars have no po- 
larization measured, so they could be either HPQs or LPQs. 
Radio galaxies are non-quasar AGN. The full sample is listed 
in Table 1. Columns 1 and 2 give the name and J2000 name for 
the source, and for some sources an alternative name is given in 
Col. 3. The coordinates of the sources are given in Cols. 4 and 
5. The start dates of the Planck scans are listed in Cols. 6 and 7. 
The average 37 GHz flux density from Metsahovi observations 
is given in Col. 8. For Col. 9 onwards, see Sect. 4. 

3. Multifrequency data 

The core of our data set is the Planck ERCSC. The con- 
struction and contents of the catalogue are described in 
Planck Collaboration (2011c). For most sources, the ERCSC 
flux density values are averages of two scans, separated by about 
six months. To enable extensive multi frequency studies with a 
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simultaneous data set, the Planck Extragalactic Working Group 
has coordinated a programme in which ground-based and space- 
home telescopes observe the sources in unison with Planck . In 
this paper we present SEDs based on the averaged ERCSC flux 
densities and on multifrequency supporting observations taken 
within two weeks of the Planck scans. In a later paper we will 
present single-epoch SEDs, constructed from all the available 
Planck and supplementary simultaneous multifrequency data. 
This collaborative programme includes 12 observatories around 
the world (Table 2). Archival data have been obtained from 
the literature and from the the search tool at the ASI (Agenzia 
Spaziale Italiana) Science Data Center (ASDC) web page 2 . 

3.1 Radio and submiliimetre data 

Centimetre-band observations were obtained with the University 
of Michigan Radio Astronomy Observatory’s (UMRAO) 26- 
m prime focus paraboloid equipped with radiometers operating 
at central frequencies of 4.8, 8.0, and 14.5 GHz. Observations 
at all three frequencies employed rotating polarimeter systems 
permitting both total flux density and linear polarization to be 
measured. A typical measurement consisted of 8-16 individ- 
ual measurements over a 20-40 min period. Frequent drift scans 
were made across stronger sources to verify the telescope point- 
ing correction curves, and observations of programme sources 
were intermixed with observations of a grid of calibrator sources 
to correct for temporal changes in the antenna aperture effi- 
ciency. The flux scale was based on observations of Cassiopeia A 
(Baars et al. 1977). Details of the calibration and analysis tech- 
niques are described in Aller et al. (1985). 

Six-frequency broadband radio spectra were obtained with 
the RATAN-600 radio telescope in transit mode by observing si- 
multaneously at 1.1, 2.3, 4.8, 7.7, 11.2, and 21.7 GHz (Parijskij 
1993; Berlin & Friedman 1996). The parameters of the receivers 
are listed in Table 3, where v c is the central frequency, Av is the 
bandwidth, AT is the sensitivity of the radiometer over a 1 s in- 
tegration, r P hys is the physical temperature of the front-end am- 
plifier, and Tgys is the noise temperature of the whole system at 
the given frequency. Data were reduced using the RATAN stan- 
dard software FADPS (Flexible Astronomical Data Processing 
System) reduction package (Verkhodanov 1997). The flux den- 
sity measurement procedure at RATAN-600 is described by 
Aliakberov et al. (1985). 

The 37 GHz observations were made with the 13.7-m 
Metsahovi radio telescope using a 1 GHz bandwidth, dual-beam 
receiver centred at 36.8 GHz. The observations were ON-ON 
observations, alternating the source and the sky in each feed 
horn. The integration time used to obtain each flux density data 
point typically ranged from 1200 to 1400 s. The detection limit 
of the telescope at 37 GHz is ~ 0.2 Jy under optimal conditions. 
Data points with a signal-to-noise ratio less than four were han- 
dled as non-detections. The flux density scale was set by obser- 
vations of DR 21. Sources NGC7027, 3C274, and 3C84 were 
used as secondary calibrators. A detailed description of the data 
reduction and analysis is given in Terasrantaet al. (1998). The 
error estimate in the flux density includes the contribution from 
the measurement rms and the uncertainty of the absolute calibra- 
tion. 

Quasi-simultaneous centimetre/millimetre radio spectra for 
a large number of Planck blazars have been obtained within 
the framework of a Fermi-GST related monitoring pro- 
gramme of gamma-ray blazars (the F -GAMMA programme, 

http://www.asdc.asi.it 


Fuhrmann et al. 2007; Angelakis et al. 2008). The frequency 
range spans 2.64 GHz to 142 GHz using the Effelsberg 100- 
m and IRAM 30-m telescopes. The Effelsberg measurements 
were conducted with the secondary focus heterodyne receivers at 
2.64, 4.85, 8.35, 10.45, 14.60, 23.05, 32.00, and 43.00 GHz. The 
observations were performed quasi-simultaneously with cross- 
scans, i.e., slewing over the source position in the azimuth and 
elevation directions, with an adaptive number of sub-scans for 
reaching the desired sensitivity (for details, see Fuhrmann et al. 
2008; Angelakis et al. 2008). Pointing offset correction, gain 
correction, atmospheric opacity correction, and sensitivity cor- 
rection have been applied to the data. 

The Institut de Radioastronomie Millimetrique (IRAM) ob- 
servations were carried out with calibrated cross-scans using the 
EMIR horizontal and vertical polarisation receivers operating at 
86.2 and 142.3 GHz. The opacity-corrected intensities were con- 
verted to the standard temperature scale and corrected for small 
remaining pointing offsets and systematic gain -elevation effects. 
The conversion to the standard flux-density scale was done us- 
ing the instantaneous conversion factors derived from frequently 
observed primary calibrators (Mars, Uranus) and secondary cal- 
ibrators (W3(OH), K3-50A, NGC 7027). From this programme, 
radio spectra measured quasi-simultaneously with the Planck 
observations have been collected for a total of 37 Planck blazars 
during the period August 2009 to June 2010, 

Many of the sources in the sample were monitored at 
15 GHz using the 40-m telescope of the Owens Valley Radio 
Observatory (OVRO) as part of a larger monitoring programme 
(Richards et al. 2011). The flux density of each source was 
measured approximately twice weekly, with occasional gaps 
due to poor weather or instrumental problems. The telescope 
is equipped with a cooled receiver installed at prime focus, 
with two symmetric off-axis corrugated hom feeds sensitive to 
left circular polarization. The telescope and receiver combina- 
tion produces a pair of approximately Gaussian beams (157" 
FWHM), separated in azimuth by 12/95. The receiver has a cen- 
tre frequency of 15.0GHz, a 3.0GHz bandwidth, and a noise- 
equivalent reception bandwidth of 2.5 GHz. Measurements were 
made using a Dicke-switched dual-beam system, with a second 
level of switching in azimuth to alternate source and sky in each 
of the two horns. Calibration was referred to 3C 286, for which a 
flux density of 3.44 Jy at 15 GHz is assumed (Baars et al. 1977). 
Details of the observations, calibration, and analysis are given 
by Richards et al. (201 1). 

The Very Large Array (VLA) and (since spring 2010) the 
Expanded VLA also observed a subset of the sources as simulta- 
neously as possible. Most of the VLA and EVLA runs were brief 
1-2 hour blocks of time. In a one-hour block of time, in addition 
to flux calibrators and phase calibrators, typically 5-8 Planck 
sources were observed. In many cases, VLA flux density and 
phase calibrators were themselves of interest, since they were 
bright enough to be detected by Planck. For these bright sources, 
the integration times could be extremely short. At 4.86 GHz and 
8.46 GHz, each target was observed for approximately 30 s. The 
integration times were typically 100 s at 22.46 GHz and 120 s 
at 43.34 GHz. All VLA/EVLA flux density measurements were 
calibrated using standard values for one or both of the primary 
calibrator sources used by NRAO (3C 48 or 3C 286), and the uv 
data were flagged, calibrated and imaged using standard NRAO 
software, AIPS or CASA. The VLA and EVLA were in different 
configurations at different times in the several months duration 
of the observations. As a consequence, the angular resolution 
changed. In addition, for a given configuration, the resolution 
was much finer at higher frequencies. For that reason, sources 
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Table 2, Optical and radio observatories that participated in the Planck multifrequency campaigns 


Radio 




Observatory 


Frequencies [GHz] 

APEX, Chile 


345 


ATCA, Australia 


4.5-40 


Effelsberg, Germany 


2.64-43 


IRAM Pico Veleta, Spain 

86, 142 


Medicina, Italy 


5, 8.3 


Metsahovi, Finland 


37 


OVRO, USA 


15 


RATAN-600, Russia 

1.1, 2.3, 4.8, 7.7,11.2,21.7 

UMRAO 


4.8, 8.0, 14.5 


VLA, USA 


5, 8, 22, 43 


Optical 

Observatory 


Band 


KVA, Spain 


R 


Xinglong, China 


i 


-600 receivers. 

f c [GHz] A/ [GHz] 

A T [mKJ r phys [K] 

T v [K] 

21.7 2.5 

3X 

15 

77 

11.2 1.4 

3 

15 

65 

7.7 1.0 

3 

15 

62 

4.8 0.9 

2.2 

15 

39 

2.3 0.4 

8 

310 

95 

1.1 0.12 

15 

310 

105 


that showed signs of resolution in any VLA configuration or at 
any VLA frequency have been carefully flagged. 

The Planck-AYCA Co-eval Observations (PACO) project 
(Massardi etal. 2011a) consists of observations of a compila- 
tion of sources selected in the Australia Telescope 20 GHz sur- 
vey (AT20G; Massardi et al. 201 lb) taken with the Australia 
Telescope Compact Array (ATCA) in the frequency range 4.5- 
40 GHz. The observations were carried out at several epochs 
close in time with Planck satellite observations covering July 
2009 - August 2010. The PACO sample is a complete, flux- 
density limited, and spectrally-selected sample over the whole 
southern sky, with the exception of the region with Galactic lat- 
itude |6| < 5°. 

The Simultaneous Medicina Planck Experiment (SiMPlE; 
Procopio et al. 2011) uses the 32-m Medicina single-dish an- 
tenna at 5 and 8.3 GHz to observe the 263 sources of the 
New’ Extragalactic WMAP Point Source (NEWPS) sample 
(Massardi et al. 2009) with 5 > 0°, and partially overlapping 
with the PACO observations for -10° < 6 < 0°. The project be- 
gan in June 2010. Because of the lack of simultaneity with the 
ERCSC, the data have been used here to match the PACO ob- 
servations in the overlapping region and to add information for 
5 > 70°, a region so far poorly covered at 5 GHz. 

Twelve sources from our sample were observed in the 
submillimetre domain with the 12-m Atacama Pathfinder 
Experiment (APEX) in Chile. The observations were made us- 
ing the LABOCA bolometer array centred at 345 GHz. Data 
were taken at two epochs in 2009: September 3-4 2009 and 
November 12 2009. The data were reduced using the script pack- 
age minicrush 3 , version 30-0ct-2009, with Uranus used as cali- 
brator. 


3 http://www.submm.caltech. edu/~sharc/cnish/ 


3.2. Optical data 

The optical observations were made with the 35-cm KVA 
(Kungliga Vetenskapsakademien) telescope on La Palma, 
Canary islands. All observations were made through the /?-band 
filter (yieff = 640 nm) using a Santa Barbara ST-8 CCD camera 
with a gain factor of 2.3 e'/ADU and readout noise of 14 elec- 
trons. The binning of pixels by 2 x 2 pixels resulted in a plate 
scale of 0798 per pixel. We obtained three to six exposures of 
1 80 s per target. The images were reduced in the standard way 
of subtracting the bias and dark frames and dividing by twilight 
flat-fields. 

The flux densities of the target and three to five stars in the 
target field were measured with aperture photometry and the 
magnitude difference between the target and a primary reference 
star in the same field was determined. Using differential mode 
makes the observations insensitive to variations in atmospheric 
transparency and accurate measurements can be obtained even 
in partially cloudy conditions. The R - band magnitude of the 
primary reference star was determined from observations made 
on photometric nights, using comparison stars in known blazar 
fields as calibrators (Fiorucci & Tosti 1996; Fiorucci et al. 1998; 
Raiteri et al. 1998; Villata et al. 1998; Nilsson et al. 2007) and 
taking into account the color term of the /?-band filter employed. 
After the /?- band magnitude of the primary reference star was de- 
termined, the object magnitudes were computed from the mag- 
nitude differences. At this phase we assumed V-R = 0.5 for the 
targets. Several stars in the field were used to check the quality 
of the photometry and stability of the primary reference. The un- 
certainties in the magnitudes include the contribution from both 
measurement and calibration errors. 

The monitoring at Xinglong Station, National Astronomical 
Observatories of China, was performed with a 60/90-cm f/3 
Schmidt telescope. The telescope is equipped with a 4096x4096 
E2V CCD, which has a pixel size of \2pm and a spatial reso- 
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lution of 173 per pixel. The observations were made with an i 
filter. Its central wavelength and passband width are 668.5 and 
51.4 nm, respectively. The exposure times are mostly 120 s but 
can range from 60 to 1 80 s, depending on weather and moon 
phase. 

3.3. Swift 

The Swift Gamma-Ray-Burst (GRB) Explorer (Gehrels et al. 

2004) is a multi-frequency space observatory devoted to the dis- 
covery and rapid follow up of GRBs. There are three instruments 
on board the spacecraft: die UV and Optical Telescope (UVOT, 
Roming et al. 2005), the X-Ray Telescope (XRT, Burrows et al. 

2005) sensitive in the 0.3 - lO.OkeV band, and the Burst Alert 
Telescope (BAT, Barthelmy et al. 2005) sensitive in the 15 - 
150keV band. Although the primary scientific goal of the satel- 
lite is the observation of GRBs, the wide frequency coverage is 
useful for AGN studies because it covers the region where the 
transition between synchrotron and inverse Compton emission 
usually occurs. 

When not engaged in GRB observations, Swift is available 
for target of opportunity (ToO) requests. The Swift team de- 
cided to devote an average of three ToO per week to the si- 
multaneous observations of Planck AGN. The simultaneity of 
Swift observations within two weeks to the Planck first or sec- 
ond scan, or both, is shown in Table 4, marked sim_lst, sim_2nd, 
or sim_lst_2nd, respectively. 

3.3.1. UVOT 

SVw’tf UVOT observations were carried out using the “filter of 
the day”, i.e., one of the lenticular filters ( V , R, U , UVW \ , 
UVXfl, and UVW2 ), unless otherwise specified in the ToO re- 
quest. Therefore images are not always available for all filters. 

The photometry analysis of all our sources was performed 
using the standard UVOT software distributed within the 
HEAsoft 6. 8. 0 package, and the calibration included in the latest 
release of the “Calibration Database” (CALDB). A specific pro- 
cedure has been developed to process all the ToO observations 
requested by the Planck project. Counts were extracted from an 
aperture of 5" radius for all filters and converted to flux densi- 
ties using the standard zero points (Poole et al. 2008). The flux 
densities were then de-reddened using the appropriate values of 
E(B - V) for each source taken from Schlegel et aL (1998), with 
A*IE(B - V) ratios calculated for UVOT filters using the mean 
galactic interstellar extinction curve from Fitzpatrick (1999). No 
variability was detected within single exposures in any filter. 

3.3.2. XRT 

The XRT is usually operated in “Auto State” mode which auto- 
matically adjusts the CCD read-out mode to the source bright- 
ness, in an attempt to avoid pile-up (see Burrows et al. 2005; 
Hill et al. 2004 for details of the XRT observing modes). As a 
consequence, part of the data sample was collected using the 
most sensitive Photon Counting (PC) mode while Windowed 
Timing (WT) mode was used for bright sources. 

The XRT data were processed with the XRTDAS software 
package (v. 2.5.1, Capalbi et al. 2005) developed at the ASDC 
and distributed by the NASA High Energy Astrophysics Archive 
Research Center (HEASARC) within the HEASoft package (v. 
6.9). Event files were calibrated and cleaned with standard filter- 
ing criteria with the xrtpipeline task, using the latest calibration 


files available in the Swift CALDB. Events in the energy range 
0.3-10keV with grades 0-12 (PC mode) and 0-2 (WT mode) 
were used for the analysis. 

Events for the spectral analysis were selected within a cir- 
cle of a 20-pixel (~ 47") radius, centred on the source posi- 
tion, which encloses about 90% of the point spread function 
(PSF) at 1.5keV (Moretti et al. 2005). For PC mode data, when 
the source count rate is above ~ 0.5 counts/s data are signifi- 
cantly affected by pile-up in the inner part of the PSF. For such 
cases, after comparing the observed PSF profile with the analyt- 
ical model derived by Moretti et al. (2005), we removed pile-up 
effects by excluding events detected within a certain inner radius 
around the source position, and used an outer radius of 30 pix- 
els. The value of the inner radius was evaluated individually for 
each observation affected by pile-up, depending on the observed 
source count rate. 

Ancillary response files were generated with the xrtmkarf 
task applying corrections for the PSF losses and CCD defects. 
Source spectra were binned to ensure a minimum of 20 counts 
per bin to allow the x 2 minimization fitting technique to be 
used. We fitted the spectra adopting an absorbed power-law 
model with photon index T x . When deviations from a single 
power-law model were found, we adopted a log-parabolic law 
of the form F(E) = (Massaro et al. 2004) which 

has been shown to fit the X-ray spectrum of blazars well (e.g., 
Giommi et al. 2005; Tramacere et al. 2009). This spectral model 
is described by only two parameters: a, the photon index at 
1 keV, and b, the curvature of the parabola. For both models 
the amount of hydrogen-equivalent column density (Afe) was 
fixed to the Galactic value along the line of sight (Kalberla et al. 
2005). 

3.4. Fermi-LAT Observations and Data Analysis 

The Large Area Telescope (LAT) onboard Fermi is an electron- 
positron pair-conversion telescope sensitive to gamma-rays of 
energies from 20MeV to above 300 GeV. The Fermi-LAT con- 
sists of a high-resolution silicon microstrip tracker, a Csl ho- 
doscopic electromagnetic calorimeter, and an anticoincidence 
detector for charged particles background identification. A full 
description of the instrument and its performance can be found 
in Atwood et al. (2009). The large field of view (-2.4 sr) allows 
the LAT to observe the full sky in survey mode every 3 hours. 
The LAT PSF strongly depends on both the energy and the con- 
version point in the tracker, but less on the incidence angle. For 
1 GeV normal incidence conversions in the upper section of the 
tracker, the PSF 68% containment radius is 0.6° . 

The Fermi-LAT data considered for this analysis cover the 
period from 4 August 2008 to 4 November 2010. They have been 
analyzed using the standard Fermi-LAI ScienceTools software 
package 4 (version v9rl6) and selecting for each source only 
photons above lOOMeV, belonging to the diffuse class (Pass6 
V3 IRF, Atwood et al. 2009) which have the lowest background 
contamination. For each source, we selected only photons within 
a 15° region of interest (Rol) centred around the source itself. 
In order to avoid background contamination from the bright 
Earth limb, time intervals where the Earth entered the LAT Field 
of View (FoV) were excluded from the data sample. In addi- 
tion, we have excluded observations in which the source was 
viewed at zenith angles larger than 105°, where Earth’s albedo 
gamma-rays increase the background contamination. The data 
were analyzed with a binned maximum likelihood technique 

4 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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(Mattox et al. 1996) using the analysis software {gtlike ) devel- 
oped by the LAT team 5 6 . A model accounting for the diffuse 
emission as well as the nearby gamma-ray sources is included 
in the fit. 

The diffuse foreground, including Galactic interstellar emis- 
sion, extragalactic gamma-ray emission, and residual cos- 
mic ray (CR) background, has been modelled using the 
models^ glLiem-v02 for the Galactic diffuse emission and 
isotropic Jem jv02 for the extra-galactic isotropic emission. Each 
source has been fit with a power law function 


dN 

dE 


N(y+\)EF 


pV+l _ F y + 1 
^max c m j n 


( 1 ) 


where both the normalization factor N and the photon index y 
are left free in the model fit. The model also includes all the 
sources within a 20° Rol included in the Fermi-LAT one year 
catalogue (Abdo et al. 2010b), modelled using power law func- 
tions. If a source included in the model is a pulsar belonging to 
the Fermi pulsar catalog (Abdo et al. 2010d), we have modelled 
the source using a power-law with exponential cut-off and the 
spectral parameters from the pulsar catalogue. For the evaluation 
of the gamma-ray SEDs, the whole energy range from 100 MeV 
to 300 GeV is divided into two equal logarithmically spaced bins 
per decade. In each energy bin the standard gtlike binned anal- 
ysis has been applied assuming for all the point sources in the 
model a power law spectrum with photon index fixed to -2.0. 
Assuming that in each energy bin the spectral shape can be ap- 
proximated by a power law, the flux density of the source in all 
selected energy bins was evaluated, requiring in each energy bin 
a test statistic (TS) 7 greater than 1 0. If the TS is lower than 1 0, an 
upper limit is evaluated in that energy bin. Only statistical errors 
for the fit parameters are shown in the plots. Systematic effects 
are mainly based on uncertainties in the LAT effective area de- 
rived from the on-orbit estimations, and are <5% near 1 GeV, 
10% below 0.10 GeV, and 20% above 10 GeV. 

The LAT gamma-ray spectra of all AGN detected by Fermi 
are studied in Abdo et al. (2010c), based on 1 1 months of Fermi- 
LAT data. For this paper, we derived the gamma-ray SEDs 
of the 104 sources in the sample in three time intervals, pre- 
sented in Table 4. In the first interval the Fermi-LAI observa- 
tions are simultaneous to the Planck first or second scan, or both, 
within two weeks (marked sim_lst, sim_2nd, or sim_lst_2nd in 
Table 4). The Fermi- LAT data have in this case been integrated 
over two weeks. In the second interval the gamma-ray data are 
quasi-simultaneous to the Planck first or second scan, or both, 
having been integrated over two months (2M_lst, 2M_2nd, or 
2M_lst_2nd in Table 4). In the third interval are sources for 
which Fermi-LAT data have been averaged over 27 months due 
to their faintness (27m in Table 4). 

Note that eight sources (0306+102, 0355+508, 0804+499, 
0945+408,11130+3815, 1413+135, 1928+738, 2005+403) that 
appear as significant gamma-ray detections in the 27 month 
dataset presented in this work, have not been included in any 
Fermi catalogue published so far (Abdo et al. 2010b). They can 
be considered new gamma-ray emitting sources, taking into ac- 


5 http://fermi.gsfc.nasa.gov/ssc/data/analysis/docunientation/Cicerone/ 
CiceroneXike lihood 

6 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ 
BackgToundModels.html 

7 The test statistic (TS) is defined as TS = -2 ln(Lo/Li) with Lq the 
likelihood of the null-hypothesis model as compared to the likelihood 
of a competitive model, L\ 


count the possibility that their association with the radio source 
is purely spatial. 

4. Spectral energy distributions 

The SEDs for the whole sample are shown in Figs. 18-121 (in 
panels on the left). The Planck measurements and the simulta- 
neous auxiliary data are indicated by red circles. The grey points 
represent archival data obtained from the literature and using the 
ASDC search tool. The measurements were regarded as simul- 
taneous if they were taken within two weeks (for radio frequen- 
cies) or five days (for optical and higher energies) of the Planck 
measurement. Previous studies of the radio variability of blazars 
(Hovatta et al. 2008; Nieppola et al. 2009) have shown that the 
two-week simultaneity limit is appropriate, as large-scale radio 
flux-density changes on that timescale are quite rare. The SEDs 
were fitted with second and third degree polynomials, namely 

log vF y = ci Gog v) 2 + c 2 0og v) + c 3 , (2) 

and 

log vF y = c 4 ( log v) 3 + c 5 Gog v) 2 + c 6 (log v) + c 7 , (3) 

where c f (i = 1 ... 7) are fit parameters. We calculated the peak 
frequencies of the components from the third degree fit by de- 
fault and used the second degree fit only if it seemed by visual in- 
spection to be more accurate. Typically, in these cases, the third 
degree fit would not yield a sensible result at all. In the SEDs in 
Figs. 18-121, third degree fits are marked with a dashed line, 
second degree fits with a dotted line. 

The fits were divided into two classes, A and B, according 
to the fit quality, A being superior. We emphasize that this divi- 
sion is the subjective view of the authors, and is only meant to 
illustrate the varying quality of the SED fits. The classes have 
not been used in calculating the correlation coefficients (see be- 
low). We were able to obtain meaningful synchrotron fits for 60 
sources (15 class A fits), and IC fits for 30 sources (10 class 
A fits). For 21 sources we were able to fit both components. 
The synchrotron and IC peak frequencies are listed in Table 1 . 
Columns 9 and 10 give the peak frequency and the fit quality for 
the synchrotron component, and Cols. 1 1 and 12 give the peak 
frequency and the fit quality for the inverse Compton compo- 
nent. 

Figure 1 shows the distributions of the logarithms of the peak 
frequencies. The synchrotron peak frequencies, logvjj, are typi- 
cally very low, as we are dealing with bright radio sources. They 
range over two orders of magnitude, from 12,2 to 14.3, the av- 
erage of the distribution being 13.2 ± 0.5. The source with the 
highest in the near-ultraviolet domain is 0716+714. The IC 
peak frequencies, log vjf, range over three orders of magnitude, 
from MeV to the high GeV gamma-ray region. The distribution 
average is 22.0 ± 0.9. The source with the highest IC peak fre- 
quency is 1 156+295, with log vjf = 23.5. The source with the 
lowest IC peak frequency, with log v* c = 20.5, is 0836+710. 

We have plotted the interdependence of the component peaks 
in Fig. 2. The two peak frequencies seem to have no significant 
correlation, which is confirmed by the Spearman rank correla- 
tion test (p = 0.301, P = 0.092, both class fits included). There 
is, however, a tendency for the sources with high vjj to have a 
high if. In the synchrotron-self-Compton (SSC) scenario, the 
separation of the component peaks depends on many factors, 
such as the electron Lorentz factors, magnetic field strengths, 
and particle densities. Therefore the lack of correlation as such 
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Fig. 1. The distributions of the synchrotron (top panel) and IC 
(bottom panel) peak frequencies. 


does not allow any strong statements to be made about the origin 
of the IC radiation (SSC, or external Compton, EC). 

We also compared our peak frequencies with those of 
Abdo et al. (2010a). These comparisons are illustrated in Fig. 3, 
where the one-to-one correspondence is shown with a dashed 
line. Both correlations are significant according to the Spearman 
test (p = 0.629, P = 0.003 for the synchrotron peaks and 
p = 0.660, P = 0.010 for the IC ones). 

5. Radio spectra 

In addition to the spectral energy distributions, we have plotted 
the standard radio spectra for the sample sources. The aim was 
to look for signs of multiple components contributing to the syn- 
chrotron emission. The spectra are shown in Figs. 18-121 (in 
panels on the right). Red circles indicate low frequency (LF; < 
70 GHz) data simultaneous to Planck , red stars indicate ERCSC 
LFI data, blue circles indicate high frequency (HF; > 70 GHz) 
data simultaneous to Planck , and blue stars indicate ERCSC HFI 
data. As can be seen, the Planck data and the supporting multi- 
frequency data provide an excellent opportunity to find the sig- 
natures of possible multiple components or anomalous spectral 
shapes. A statistical study of extragalactic radio sources in the 
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Fig. 2. The correlation of the peak frequencies of the synchrotron 
and the IC components. Solid circles denote class A fits, and 
open circles class B fits. 
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Fig. 3. Comparison between our synchrotron and IC peak fre- 
quencies and those of Abdo et al. (2010a). The synchrotron peak 
frequencies are shown in the top panel and IC peak frequencies 
in the bottom panel. The dashed lines illustrate the one-to-one 
correspondence. Solid circles denote class A fits, and open cir- 
cles class B fits. 
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ERCSC is presented in Planck Collaboration (201 li), and SEDs 
of extreme radio sources, such as Gigahertz-Peaked Spectrum 
(GPS) sources in Planck Collaboration (201 lj). 

In Figs. 4 and 5 we have plotted the distributions of the LF 
and HF spectral indices using only the Planck data for the latter. 
When plotting the HF a values, we have only taken into account 
those sources with three or more data points. This leaves us with 
84 sources out of the 104. The sources with only one Planck 
scar are shown hatched in Figs. 4 and 5. Figure 6 shows otlf 
versus o^f- As expected, the LF indices are fairly flat, with an 
average of -0.06. Their distribution is narrow, with 91% of the 
indices being in the range a L F = -0.5-0.5. There are a couple 
of sources with remarkably steep LF spectra, namely 0552+398 
and 2021+614, both having a L f < -0.8. At the other end we 
have 0007+106 which has an inverted spectrum with (*lf = 0.86. 
Also 1 228+1 26 has a steep spectrum, although the fit value may 
be exaggerated in this case (see Fig. 75). For clarity, this source 
has been omitted from Figs. 4 and 6. 

The HF indices are concentrated around -0.5, having an av- 
erage of -0.56. Rather similar distributions for LF and HF spec- 
tral indices were found in Planck Collaboration (201 li), espe- 
cially in their Fig. 7. However, a large part of the othf distribu- 
tion is in the flat domain with qthf > -0.5. Extreme examples are 
1413+135 (<*h f = 0.02) and 0954+658 (ot H f = 0.34). It is unex- 
pected that the spectral indices of blazars at sub-millimetre and 
FIR frequencies should be this flat. Only 28 of the 84 sources 
with HF spectral fits shown in Figs. 18-121 have a < -0.7, 
the ‘canonical” value for optically thin spectra, corresponding 
to an electron energy index s of about 2.4. This has two possible 
explanations: either the total HF spectra are defined by several 
underlying components, or the energy spectrum of the electron 
population producing the radiation is much harder than generally 
assumed. Although some HF spectra show evidence for several 
sub- components or even an upturn at the highest frequencies, 
many others appear to be perfectly straight and therefore per- 
haps originate from a single optically thin component. At the 
highest Planck frequencies the lifetimes of radiating electrons 
are short, and one would in most cases expect to see the opti- 
cally thin spectra steepened by energy losses, with indices rang- 
ing from the original orthin to orthin - 0.5 and beyond. This makes 
spectral indices > -0.7 even more remarkable. 

To rule out contamination by dust in the host galaxies, we 
have estimated the possible contribution of an extremely bright 
submillimetre galaxy. We also estimated the flux densities of 
2251+158 at several submillimetre bands in the case its lumi- 
nosity were similar to the submillimetre galaxy. The compari- 
son shows that if the dust component in 2251 + 158 were of the 
same brightness as in the extremely bright submillimetre galaxy, 
its emission would be completely swamped by the non-thermal 
emission from the jet. The LF and HF spectral indices have 
been calculated using ERCSC data which for most sources is 
an average of two Planck scans separated by about six months. 
However, the HF spectral index distribution for the sources that 
have been observed only once by Planck is quite similar to that 
of the whole sample. The implications of this result are discussed 
in more detail in Sect. 7. 

As Fig. 6 shows, BLOs, HPQs and LPQs have different o?hf 
distributions. In Table 5 we have listed the average spectral in- 
dices for both frequency ranges for all subgroups. The LF spec- 
tra of radio galaxies are very steep on average, but the value 
is greatly influenced by the uncertain fit of 1228+126. If this 
source is left out of the calculations, the average otlf for galax- 
ies is 0.10 ± 0.52. LPQs have the steepest spectra in the HF 
range. However, the standard deviations of the samples are sub- 


Table 5. Average spectral indices for AGN subclasses. The er- 
rors correspond to one standard deviation. 


Class 

number LF 

number HF 

average ar LF 

average « H f 

BLO 

21 

20 

-0.01 ± 0.25 

-0.34 ± 0.26 

HPQ 

33 

28 

-0.06 ± 0.22 

-0.58 ±0.27 

LPQ 

21 

18 

-0.08 ± 0.30 

-0.80 ±0.27 

QSO 

19 

10 

-0.05 ± 0.28 

-0.57 ± 0.26 

GAL 

9 

7 

-0.25 ± 0.66 

-0.55 ±0.19 



Fig. 4. Distribution of LF spectral indices for the whole sample 
(1 04 sources). The sources that have only been scanned once by 
Planck are shown hatched (46 sources). 

stantial, which can also be seen as the large scatter in Fig. 6. 
To quantify the possible differences between the AGN classes, 
we ran the Kruskal- Wallis ANOVA test. For the LF indices we 
could find no significant differences. However, in the case of the 
HF indices, the distributions of the subgroups differ significantly 
(P < 0.001). Multiple-comparisons test with t distribution tells 
us that LPQs differ from all other subgroups (P < 0.031). Also, 
BLO indices have a different distribution from all quasar sub- 
groups (P < 0.019). 

6. Modelling the SEDs and spectra of blazars 

Spectral energy distributions rapidly became one of the main 
tools for understanding blazar physics after the EGRET instru- 
ment aboard the Compton Gamma-ray Observatory (CGRO) 
satellite discovered strong gamma-ray emission from 3C279 
(Hartman et al. 1992) and, subsequently, from a large number of 
other radio-bright AGN. It was rapidly accepted that the double- 
peaked overall shape of the radio-to-gamma SED was due to 
synchrotron and inverse Compton radiation, but beyond this the 
agreement ended. The data simply lacked the sensitivity and 
the time coverage to sufficiently constrain the models. Even the 
most detailed and convincing effort for 3C279 (Hartman et al. 
2001) could be criticized for unrealistic physical assumptions 
(Lindfors et al. 2005). 

With the coming of the Fermi and Swift satellites, ground- 
based TeV telescopes, and other satellite data, the inverse 
Compton high-energy SED of a large number of AGN can now 
be determined with hugely improved accuracy and time cov- 
erage. This has initiated a new era of AGN modelling (e.g., 


8 



Planck Collaboration: SEDs and radio spectra of northern AGN 



<*hf 


Fig. 5. Distribution of HF spectral indices for the 84 sources that 
had three or more data points. The sources that have only been 
scanned once by Planck are shown hatched (40 sources). 



Fig. 6. LF spectral indices qtlf versus HF spectral indices c*hf 
for the various AGN subclasses. 


Proceedings of the Workshop “Fermi meets Jansky: AGN in 
Radio and Gamma Rays”; Fermi meets Jansky 2010). However, 
since the IC spectrum is created by the relativistic electrons scat- 
tering either the synchrotron photons (SSC) or external photons 
(EC), accurate knowledge of the shape of the synchrotron com- 
ponent is essential for all realistic modelling. Here the Planck 
data are invaluable, since they provide the first simultaneous 
multifrequency coverage of the crucial radio-to-IR part of the 
synchrotron spectrum. 

However, most of the ERCSC data arc not suitable for de- 
tailed modelling, since they are averages of two spectra taken 
about six months apart. Considering the strong SED variabil- 
ity, accurate modelling can be done only with the final Planck 
datasets and their several individual scans. Here we limit our- 
selves to a short discussion of how Planck and supporting data 
can be used to improve our understanding of blazar physics. 

The basic picture of a blazar is well established: a black 
hole surrounded by an accretion disk and a broad-line emit- 


ting region (BLR), and a relativistic jet which produces the 
lower-frequency bump of the SED through synchrotron emis- 
sion. The jet is not totally stable; VLBI imaging shows how 
new components, identified as relativistic shocks in the jet, 
emerge from the radio core at intervals ranging from months 
to years. These growing and decaying shocks are also respon- 
sible for the total flux density variations in the radio-to-IR, 
and at least partly also in the optical regime. The important 
thing to note is that the synchrotron SED is rarely, if ever, the 
product of a single synchrotron-radiating component. The na- 
ture of the radio core is still being debated, as is the nature 
and the extent of the jet flow upstream of the radio core, be- 
tween it and the black hole (e.g., Marscher & Jorstad 2010; 
Marscher et al. 2010). Shocked jet models (Marscher & Gear 
1985; Hughes etal. 1989; Valtaoja et al. 1992) provide a good 
account of the multifrequency variations in the total flux density. 
In the particular case of 3C279, Lindfors et al. (2006) showed 
that the optical variations can also be moderately well explained 
as the higher-frequency part of the shock emission. 

The origin of the IC component is a matter of intense de- 
bate. In the CGRO era, the favoured alternative was that the 
high-frequency radiation originates close to the black hole and 
the accretion disk, inside the BLR and well upstream of the ra- 
dio core. With a strong radiation field external to the jet, exter- 
nal Compton is then likely to be the dominant emission mech- 
anism. A minority view was that at least a part of the gamma- 
ray flaring comes from much farther out, at or downstream 
from the radio core. The evidence for this was that, statisti- 
cally speaking, strong gamma-ray flaring tended to occur af- 
ter the beginning of a millimetre flare (Valtaoja & Terasranta 
1996; Lahteenmaki & Valtaoja 2003) and after the ejection of 
a new VLBI component (Jorstad et al. 2001). Occurring well 
outside the BLR, SSC seems in this case to be the only vi- 
able emission mechanism, but on the other hand simple SSC 
models have generally failed to explain the gamma-ray emis- 
sion, requiring unrealistically low magnetic field strengths in 
other than the extreme TeV BL Lac objects. With new Fermi 
data, the evidence for gamma-ray emission from the radio jet 
has strengthened (Kovalev et al. 2009; Pushkarev et al. 2010; 
Leon-Tavares et al. 201 1 ; Nieppola et al. 2010; Tornikoski et al. 
20 1 0; Valtaoja et al. 20 1 0), but the source of the upscattered pho- 
tons remains a puzzle. 

Most of the attempts to model the high-energy bump have 
adopted a “reverse engineering” method, starting from the sec- 
ondary IC SED and not from the primary synchrotron SED. 
The original electron energy spectrum is considered to be rather 
freely adjustable, with parameters chosen to produce a good fit 
of the calculated IC spectrum to the observations within the as- 
sumed model. In most cases, the synchrotron spectrum is as- 
sumed to originate in a single homogeneous component, and of- 
ten the spectrum is adjusted only to the optical data, without 
any attempt to explain the radio-to-IR part of the SED. There 
is not much observational or theoretical evidence for the actual 
physical existence of such an IR-to-opdcal component: it is sim- 
ply assumed to be there, because it can produce the observed IC 
spectrum. 

To some extent such an approach is justifiable if the gamma- 
rays originate upstream of the radio core, and are therefore not 
directly connected to the observed radio jet and shocks, and their 
lower-frequency synchrotron emission. However, even in this 
case physically more realistic models must take into account also 
the radio-to-IR part, for several reasons. First, whatever the exact 
mechanism, the radio shocks also contribute to the IC spectrum, 
and their contribution must be accounted for in the model fits. 
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Assuming reasonable physical parameters, the shocks can quite 
often produce part or all of the observed X-ray emission through 
the SSC mechanism (e.g., Lindfors et al. 2006). Indeed, a stan- 
dard method for estimating the physical parameters of VLBI 
components is to assume that the X-ray emission is produced 
through SSC (Ghisellini et al. 1993; Marscher 1987). Secondly, 
the components seen in the radio, especially in the Planck HFI 
regime, contribute to the optical synchrotron emission, and this 
contribution must be subtracted from the putative IR-to-optical 
component, which, as explained above, in many models gener- 
ates the whole IC spectrum. For example, in the present sam- 
ple, at least 0235+164, 0917+449, 0954+658, 2227-088, and 
especially 2251 + 158 show ongoing millimctre-to-submillimetre 
flaring, which leaves little room for any additional components 
between the highest Planck frequencies and the optical. This is 
the case in particular if the optically thin spectra are originally 
very flat, as we argue in Sects. 6.1 and 7. 

On the other hand, if gamma-rays do come from the shocks 
at, or downstream from, the radio core, as increasing evidence 
indicates, the obvious first step is to model the shocks them- 
selves and their spectra as accurately as possible using Planck 
and other radio data, and only after that attempt to calculate the 
IC radiation they produce. 

Finally, it is clear that the more realistic multicomponent 
models, whatever their exact details, will give rather different 
physical parameter values than the single-component models. As 
an example, if the X-ray emission and the gamma-ray emission 
are not produced by a single component but by two (or even 
more) components, the derived electron energy spectral parame- 
ters will be quite different. 

6. 1. Some examples of different spectra and their 
implications 

In this section we highlight five examples of different spectra of 
ERCSC sources which have been scanned only once by Planck 
and therefore have true snapshot spectra presented in Figs. 1 8- 
121. As stated before, our aim at this stage is not to model these 
sources, but to draw attention to some features commonly seen 
in the spectral energy distributions, and to their implications for 
future multi-epoch modelling of the Planck sources. We have 
fitted the radio-to-optical spectra with one or several standard 
self-absorbed sychrotron components (based on a power-law as- 
sumption for the electron energy spectrum). These are not model 
fits from a numerical, physical code (except for Figs. 16 and 17), 
but rather are meant to guide the eye and to illustrate the general 
shape of the radio-to-optical spectra. 

We also use support observations to address the behaviour 
of the five example sources. In Figs. 7-11 we have plotted the 
37 GHz long-term flux density curves for the five sources. The 
long-term data have been taken at Metsahovi Radio Observatory. 

Even without numerical modelling we can address two im- 
portant topics with these sources. The first one is the spectral 
flatness of the optically thin radiation possibly hinting at non- 
standard acceleration processes. The second one concerns the 
origin of the IC component. As Figs. 12-17 show there does not 
appear to be room for additional synchrotron components be- 
tween the Planck and the optical frequencies. The IC radiation 
must therefore originate in the same shocks that also produce the 
radio emission, in contradiction to most theoretical models. 

0234+285 (Fig. 12). This source was in a very uneventful 
state during the Planck observing period (Fig. 7) and shows a 
spectrum typical of sources in a quiescent stage. The radio spec- 
trum can be modelled with a single synchrotron component, hav- 
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Fig. 7. The 37 GHz iightcurve of 0234+285, measured at 
Metsahovi Radio Observatory, showing continuous variability. 
The period when Planck was observing the source is indicated 
by the hatched region in the lower panel, which is a blow-up of 
the last year. 
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Fig. 8. The 37 GHz Iightcurve of 0235+164, details as in Fig. 7. 


1253 055 




P rl -i- , T .J 

Li # 6 *•* * 

r xW sr * f * 

fr i w « 

; i 


' 



2009.5 2010 2010.E 

Years 
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in Fig. 7. 
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Fig. 10. The 37 GHz lightcurve of 141 3+ 135, details as in Fig. 7. 
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Fig. 11. The 37 GHz lightcurve of 2251 + 158 (3C454.3), details 
as in Fig. 7. 


ing arthin = -0.61, which apparently steepens by Aar = -0.5 
and smoothly joins the optical spectrum (sec also Fig. 29). This 
synchrotron component can be identified with the latest signifi- 
cant flare (new shock), which peaked about two years earlier at 
37 GHz (Fig. 7). During low flux levels, in the absence of recent 
flares/shocks, the underlying, relatively stable jet may also con- 
tribute significantly to the synchrotron emission. Another similar 
example is the source 0007+106 (Fig. 19). Lacking mid-IR data, 
we cannot prove that there is no additional synchrotron compo- 
nent, between the Planck frequencies and the optical, that might 
be the main source of the IC flux (Fig. 29). However, the most 
straightforward interpretation of the data is that the 1C flux also 
comes from the synchrotron component outlined in Fig. 12: from 
a ra±er old shock far downstream the jet, and far beyond the 
BLR. 

0235+164 (Fig. 13). The bumpiness of the radio spectrum 
shows clear evidence of multiple components (shocks) con- 
tributing to the observed total flux density. We have outlined two 
possible synchrotron components, which together could explain 
most of the observed radio spectrum. Although 0235+164 was 
relatively quiescent during the Planck observation (Fig. 8), there 
were at least four strong flares during the previous five years, and 
they all are likely to be still contributing significantly to the total 
radio flux density, at least around 10 10 Hz and below. Assuming 


reasonable physical parameters for these shocks, they can also 
produce significant amounts of X-rays through the SSC mecha- 
nism, which must be accounted for in a proper multicomponent 
modelling of the IC spectrum. As with 0234+285, the higher- 
turnover component with arthin = -0.4 can be joined smoothly to 
the optical spectrum assuming a steepening by about Aar = -0.5, 
leaving no room for additional mid-IR components. Fermi mea- 
sured only upper limits for the gamma-ray flux, which may be 
due to the fact that even the most recent shock, peaking around 
100 GHz, is already about two years old (Fig. 8). 

1253-055 (3C279). (Fig. 14). Again, the highest-peaking 
radio component with arthin = -0.6 can be made to smoothly 
join the steep optical spectrum. From Fig. 14 it is rather obvious 
that there can be no additional mid-IR components which might 
provide a significant contribution to the IC spectrum. If we as- 
sume that there is another synchrotron component turning over 
somewhere above the Planck frequencies and joining smoothly 
to the optical spectrum, Fig. 14 shows that its flux density can- 
not exceed the dashed line at any IR frequency between 10 12 and 
10 14 5 Hz. This is the case even if we make the assumption that 
the spectrum of the component seen at the highest Planck fre- 
quencies steepens very rapidly above 857 GHz. Consequently, 
the Fermi gamma-rays must originate in the latest shock, which 
started a few months before the Planck observation (Fig. 9). 

1413+135 (Fig. 15). This source exhibits a totally flat spec- 
trum, with both the LF and the HF spectral indices approxi- 
mately 0. Possible explanations are: an underlying electron en- 
ergy index s ~ l, which is rather unlikely; die superposition 
of many self-absorbed synchrotron components (shocks of var- 
ious ages) with different turnover frequencies; or a continuous 
jet dominating the flux and having an integrated spectral index 
close to zero over a wide frequency range. As Fig. 10 shows, the 
source has been rather inactive for over a decade, with no ma- 
jor ourbursts. The HFI spectrum may indicate an upturn, hinting 
at the presence of an IR component. We show a possible spec- 
tral decomposition, including an IR-peaking component. Further 
Planck observations are needed to ascertain the reality and the 
nature of this component, which may be either a synchrotron 
flare seen at an early developmental stage or, possibly, an in- 
frared dust component. The correct alternative can be identified 
through variability information. 

2251+158 (3C 454.3). (Figs. 16 and 17). This source ex- 
hibits the most spectacular case of ongoing IR flaring, show- 
ing a very strong synchrotron component with a self-absorption 
turnover at about 80 GHz and a very flat (or = -0.2) optically 
thin spectrum. As Fig. 1 1 shows, a series of spectacular flar- 
ing events has occurred in 2251 + 158 during the last five years, 
with Planck observing the source during the early stages of the 
strongest outburst ever seen in the source. In Fig. 1 6 we show 
a model fit to the radio-to-optical data, using a preliminary ver- 
sion of the code that we are developing (Tammi et al., in prepa- 
ration). A close-up of the fit in the radio regime is shown in 
Fig. 17. As Fig. 119 also shows, it is impossible to squeeze in 
another component between the strong radio-IR flare and the 
optical regime. The gamma-ray emission must come from the 
growing shock that we see both in the total flux density monitor- 
ing curve (Fig. 11) and in the SED (Fig. 119). Most of the theo- 
retical model scenarios assume that gamma-rays originate close 
to the black hole and the accretion disk, well within the BLR and 
far upstream of the radio core. This is not a viable alternative for 
225 1 + 1 58, in which gamma-rays must come from the new shock 
component, parsecs downstream from the radio core. As Fig. 1 1 
shows, at the time of the Planck observation the shock had al- 
ready been growing for about half a year, which, changed into 
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Fig. 12. A possible single-component spectrum for 0234+285 to 
illustrate the general shape. Red circles, LF data simultaneous to 
Planet, blue circles, HF data simultaneous to Planck . 
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Fig. 13. A possible multiple-component spectrum for 0235+164 
to illustrate the general shape. Symbols as in Fig. 12. 


the source frame, translates into a distance L = D 2 c(l +z) _1 A/ Dbs 
(where D is the Doppler boosting factor, c the speed of light, 
z the redshift, and A^bs the time elapsed from the onset of the 
flare), of at least 10 pc down the jet, assuming a modest value of 
10, or 90 pc using the value of D = 33.2 from Hovatta et al. 
(2009). This is the distance that the shock has moved down- 
stream after emerging from the radio core. 


7. The shape of the synchrotron spectrum: 
nonstandard acceleration mechanisms? 

If the electron energy distribution is assumed to be a simple 
power-law with the form N(E) = KE~ S , the synchrotron spec- 
trum of a single homogeneous source has an optically thick low- 
frequency part with a spectral index a = +2.5, a turnover, and an 
optically thin part with a spectral index or related to the electron 
energy index by a = -( s - l)/2. Energy losses eventually cause 
the spectrum to steepen to a - 0.5. 


1253-055 
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Fig. 14. The radio spectrum of 1253-055 shows that there is no 
room for another component between radio and optical regimes. 
Symbols as in Fig. 12. 
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Fig. 15. The radio spectrum of 1413+135 is very flat and can be 
fitted with multiple components. Symbols as in Fig. 12. 
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Fig. 16. The radio-optical SED of 2251+158 (3C454.3) fitted 
with one component and synchrotron emission from the jet 
(dashed line) using a model by Tammi et al. (in preparation). 
The radio spectrum is shown in Fig. 17. 
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Fig. 17. A close-up of the radio spectrum and model fit of 
2251+158 (3C 454.3), details as in Fig. 16. 


Since the earliest times of AGN research, the canonical value 
for or has been assumed to be about -0.75, corresponding to s = 
2.5. Such an electron energy index emerges naturally from the 
simplest first-order Fermi acceleration mechanism, and blazars 
also tend to show slopes around this value at the high-frequency 
(10 to 100 GHz) end of the radio spectrum, where the sources 
are typically assumed to be optically thin. However, Planck data 
now challenge this commonly accepted view. As Fig. 4 shows, 
the LF spectral slopes are, in the main, quite flat, and sometimes 
even rising. It is likely that in most cases this results from the 
superposition of several distinct synchrotron components, since 
the spectra are generally not smooth, and we also know from 
VLBI observations that there are nearly always several compo- 
nents (shocks) present in any source (see also Sect. 6). The HF 
spectra (Fig. 5), on the other hand, are often straight within the 
Planck flux density error bars, and give a definite impression 
of optically thin radiation. In some cases, such as 2251+158, 
the spectrum even shows a clear self-absorption turnover at high 
frequencies, followed by a flat optically thin part with a = -0.2 
(Figs. 16 and 17). 

The flatness of the HF spectra, shown in Fig. 5, is remark- 
able. Only a third of the sources have a steeper than -0.7, while 
15 have a > -0.3. While the present data are not sufficient to ex- 
clude the possibility that multiple components produce the flat- 
ness, and some sources such as 0954+658 do show evidence of 
an additional IR component (see also, e.g., Raiteri et al. 1999), 
we consider such an explanation unlikely for the whole sample. 
The errors in o?hf cannot explain the very flat spectral indices, 
since they generally are smaller than Aar % 0.1. 

In Table 6 we list the 10 sources that fullfill the following 
criteria. They have been scanned only once by Planck and they 
have or H F > -0.5 + Aar, where Aa is the HF spectral fit error 
shown in Figs. 18 - 121. These form the most extreme end of 
the othf distribution shown in Fig. 5. Three of them may show 
indications of multiple HF components in the spectra, but for the 
seven other sources the HF spectra appear to be straight within 
the Planck errors. If this is indeed the case, they cannot be ex- 
plained with standard acceleration mechanisms with s > 2. 

The rise of a submillimetre flare is quite fast (see, for exam- 
ple, the JCMT monitoring data of Stevens et al. 1994); the self- 
absorption turnover of a fresh shock component passes through 
the Planck HFI frequencies rather rapidly, and one therefore 
sees it only rarely in the snapshot high frequency spectra. (In 
the present data only 2251 + 158 is clearly such a case, al- 


Table 6. A list of the 10 sources scanned only once by Planck 
and with qthf > -0.5+ Aar. S = straight HF spectrum; B = bumpy 
HF spectrum, possibly indicating multiple components. 


Source 

nr hf 

Spectral shape 

0003-066 

-0.38 ±0.09 

S 

0059+581 

-0.19 ±0.09 

S 

0235+164 

-0.40 ± 0.03 

S 

0430+052 

-0.34 ± 0.05 

s 

1308+326 

-0.41 ± 0.08 

B? 

1413+135 

0.02 ±0.12 

B 

1418+546 

-0.16 ±0.04 

S 

1652+398 

-0.21 ± 0.04 

B 

1823+568 

-0.31+0.13 

S 

2251+158 

-0.20 ± 0.02 

S 


though 0235+164 and 1652+398 might be others.) Instead, in 
most cases we expect to see optically thin submillimetre spec- 
tra which have already experienced energy losses. Therefore the 
observed high-frequency spectral- index distribution of a sam- 
ple of sources should show a range from ar t hin = -(s - l)/2 to 
<*thin -0.5 and beyond, as the exponential steepening sets in, with 
spectral indices steeper than the initial arthin dominating. The ob- 
served distribution in Fig, 5 is thus incompatible with an initial 
athin around -0.7, as we should then see a few values around 
a = -0.7 and a distribution peaking towards -1.2, reflecting the 
A a = -0.5 steepening due to energy losses. 

Instead, the spectral index distribution is compatible with an 
electron index s « 1.5, resulting in a distribution of or^in hav- 
ing smallest values around -0.2 to -0.3 and a maximum around 
-0.7. Such a suggestion of a very hard original electron spec- 
trum is not entirely new; Valtaoja et al. (1988) studied the shock 
spectra and concluded that a is around -0.2, and Hughes et al. 
(1991) also concluded that the flux density variations can best 
be modelled with very hard spectra. Gear et al. (1994) also found 
that submillimetre spectra in a sample of 48 sources tended to be 
flat. Indeed, the distribution they found (their figure 7) is quite 
similar to ours (Fig. 5). 

As noted above, alternative explanations for the othf distri- 
bution cannot be totally excluded. However, we consider sources 
such as 2251+158 (Figs. 16 and 17) to be compelling evidence 
for the existence of very hard non-standard electron spectra. The 
implications of this possibility are discussed below. 

If the original a thin really is around -0.2, the distribution 
shown in Fig. 5 becomes easy to understand: in many slightly 
older HFI components we see the expected spectral steepening 
from -0.2 to -0.7 and beyond. One way to test this hypothesis is 
to try to estimate the ages of the highest- frequency spectral com- 
ponents, either from their turnover frequencies or from the total 
flux-density monitoring: the older the component, the lower the 
turnover frequency and the steeper the spectra. However, this 
analysis is beyond the scope of this paper and reliable conclu- 
sions require more Planck data than presently exist. 

An electron energy spectral index as hard as 1.5 has ma- 
jor implications for the acceleration mechanisms dominating in 
blazar shocks. Particle acceleration efficiency is strongly gov- 
erned by the compression ratio of the flow, which in turn is con- 
fined by well-known plasma physics, preventing the spectral in- 
dex from becoming arbitrarily small. For a relativistic shock the 
traditional first-order Fermi mechanism usually produces spec- 
tral indices of 2.2 or so (and not 2.5, which follows from the 
nonrelativistic theory), but not smaller. There are a few options, 
however, for bypassing the s « 2.2 limit. For example, power- 
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law spectra with s £ 1 have been found in certain kinds of 
oblique shocks (e.g., Kirk & Heavens 1989), as well as in par** 
allel shocks in the case of large-angle particle scattering due 
to very strong turbulence (e.g., Stecker et al. 2007), or when 
the compression ratio felt by the particles is higher than that 
felt by the plasma (Ellison et al. 1990; Virtanen & Vainio 2005). 
The last alternative is an especially interesting one, because 
the combination of low density and relatively strong magnetic 
field — conditions likely to be found in Poynting-flux-dominated 
AGN jets — quickly leads to non-negligible AlfVen speeds and 
enhancement of the compression ratio felt by the particles and, 
consequently, to very hard particle spectra (Tammi 2008 a, b, and 
references therein). High Alfven speeds and scattering centres 
that are not frozen into the plasma also enable second-order 
Fermi acceleration, which can produce power-law spectra even 
with s < 1 in relativistic parallel shocks (Virtanen & Vainio 
2005) and can work on time-scales comparable to the fastest 
blazar flares (Tammi & Duffy 2009). 

.As a simple order-of-magnitude reality check we can esti- 
mate the conditions required by the standard first-order accel- 
eration in a parallel step shock in the presence of an increased 
scattering-centre compression ratio due to turbulence trans- 
mission (Vainio et al. 2003). Following the analysis of Tammi 
(2008b), we can estimate that in order to achieve s « 1.5 in a 
shock with Lorentz factor F = 10, magnetic field strength of the 
order of 0.02-2 Gauss (depending on the composition and the 
density of the plasma) is needed. This is in good agreement with 
the usual modelling parameters of these sources, although we 
emphasise that the analysis used here is a simple one and may 
be limited to parallel shocks and weak turbulence. 

8. Conclusions 

We have presented the averaged Planck ERCSC spectra together 
with supporting ground and satellite observations obtained si- 
multaneously with the Planck scans. The 104 SEDs, supple- 
mented with archival data, demonstrate the usefulness of Planck 
data in determining the SED peak frequencies and fluxes and in 
modelling the spectral energy distributions in greater detail. In 
particular, the data demonstrate that the synchrotron spectrum 
contains contributions from several physically distinct AGN 
components: the LF and HF spectra are rarely smooth, except at 
the highest radio frequencies where the source component with 
the highest turnover becomes optically thin. The LF spectral in- 
dices cluster around or = 0, also indicating the superposition of 
many components with different turnover frequencies. In mod- 
elling the synchrotron and inverse-Compton SEDs, a multicom- 
ponent approach is therefore necessary. While physical mod- 
elling is beyond the scope of this paper, which uses averaged, 
rather than single-epoch, Planck ERCSC flux densities, we have 
shown some examples of multicomponent decomposition of the 
SEDs and the conclusions that can be drawn from the data. 

The sources display remarkable variability, which must also 
be taken into account in modelling the SEDs. We have presented 
some examples of total flux density monitoring to demonstrate 
how sources have been observed during Planck scans in different 
stages of activity. 

The Planck HFI spectra are remarkably flat, with only a mi- 
nority of sources having spectral indices steeper than the “canon- 
ical'' a = -0.7, and many apparently having or^n around -0.2 
to -0.4. We suggest that the most likely interpretation for the 
veiy flat and straight synchrotron spectra observed by Planck is 
a very hard original electron spectral energy index, clearly be- 
low s = 2. While this possiblity must be confirmed by fiiture 


data, an index s = 1.5 is compatible with the HF spectral index 
distribution seen in our sample. In a number of sources, energy 
losses steepen the HF spectra by Aa = 0.5 from about -0.3 to 
about -0.8, and further to -1 and beyond, which can explain 
the observed distribution. The hardness of the electron spectrum 
would also mean that the synchrotron SED peak is in most cases 
related to the maximum electron Lorentz factor achieved in the 
acceleration process. The hard spectra also require a rethinking 
of the acceleration processes dominating in the relativistic jets 
and shocks. 

For some sources the HF spectra are flat or even rising, in- 
dicating the presence of multiple synchrotron -emitting compo- 
nents peaking in the gigahertz to terahertz regime, and ongo- 
ing flaring at still higher frequencies. However, in many cases 
the Planck HFI spectrum appears to join the optical spec- 
trum smoothly, leaving little room for additional still higher- 
frequency components. In these cases, the inverse-Compton 
gamma-rays must originate in the synchrotron (i.e., shock) com- 
ponents seen in the Planck data. 
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Table 1. The complete I Jy northern sample of AGN. Column 8 refers to the 37 GHz observations of Metsahovi Radio Observatory. 
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2010-01-24 


1.64 





0212+735 

0217+7349 


02:17:30.81 

73:49:32.60 

2009-09-13 

2010-02-11 

2.05 

12.4 

B 

... 


0224+671 

0228+6721 


02:28:50.05 

67:21:03.00 

2009-09-13 

2010-02-10 

1.28 


... 

... 


0234+285 

0237+2848 

4C28.07 

02:37:52.41 

28:48:09.00 

2010-02-03 


2.86 

13.1 

A 

22,1 

B 

0235+164 

0238+1636 


02:38:38.80 

16:36:59,00 

2010-02-01 


2.33 

13.1 

A 

... 


0238-084 

0241-0815 


02:41:04.80 

-08:15:20.75 

2010-01-24 


1.33 





0306+102 

0309+1029 

PKS0306+102 

03:09:03.60 

10:29:16.00 

2010-02-06 


1.01 

12.9 

B 



0316+413 

0319+4130 

3C84 

03:19:48.16 

41:30:42.10 

2010-02-13 


14.01 

... 

... 

22.9 

B 

0333+321 

0336+3218 

NR AO 140 

03:36:30.11 

32:18:29.30 

2010-02-15 

... 

1.49 



-♦* 


0336-019 

0339-0146 

CTA026 

03:39:30.94 

-01:46:36.00 

2010-02-10 


2.21 

13.1 

B 

... 


0355+508 

0359+5057 

NRAO150 

03:59:29.75 

50:57:50.20 

2009-09-13 

2010-02-19 

5.92 


... 

... 

... 

0415+379 

0418+3801 

3C111 

04:18:21.28 

38:01:35.80 

2010-02-22 


5.97 

... 

... 

... 


0420-014 

0423-0120 

OA129 

04:23:15.80 

-01:20:33.10 

2010-02-21 

»»♦ 

7.27 

13.2 

B 



0430+052 

0433+0521 

3C120 

04:33:11.10 

05:21:15.60 

2010-02-24 


2.65 





0446+112 

0449+1121 

PKS0446+112 

04:49:07.67 

11:21:28.00 

2010-02-27 


1.41 

IM 

iM 

1|# 


0458-020 

0501-0159 

PKS0458— 020 

05:01:12.81 

-01:59:14.00 

2010-03-01 


1.86 

... 

... 

... 


0507+179 

0510+1800 


05:10:02.37 

18:00:41.58 

2009-09-13 

2010-03-03 

1.14 

|it 

l 

•it 


0528+134 

0530+1331 

PKS0528+134 

05:30:56.42 

13:31:55.15 

2009-09-14 

2010-03-07 

4.55 



21.8 

B 

0552+398 

0555+3948 

DA 193 

05:55:30.81 

39:48:49.17 

2009-09-22 

2010-03-09 

3.18 

... 

... 

21.4 

B 

0605-085 

0607-0834 

PKS0605-085 

06:07:59.70 

-08:34:50.00 

2009-09-19 

2010-03-17 

1.73 

13.2 

B 

... 


0642+449 

0646+4451 

OH471 

06:46:32.03 

44:51:16.59 

2009-09-30 

2010-03-16 

2.36 

12.5 

B 

22.7 

B 

0716+714 

0721+7120 


07:21:53.30 

71:20:36.00 

2009-10-04 

2010-03-10 

2.12 

14.3 

A 

23.4 

B 

0723-008 

0725-0054 

PKS0723-008 

07:25:50.70 

-00:54:56.00 

2009-10-09 

2010-04-02 

1.05 

13.0 

A 

... 


0735+178 

0738+1742 

PKS0735+17 

07:38:07.40 

17:42:19.00 

2009-10-10 

2010-03-31 

1.93 



... 


0736+017 

0739+0137 


07:39:18.03 

01:37:04.60 

2009-10-12 

2010-04-04 

1.85 

14.0 

B 

... 


0748+126 

0750+1231 


07:50:52.05 

12:31:04.83 

2009-10-14 

2010-04-04 

2.42 

12.7 

A 

21.5 

A 

0754+100 

0757+0956 

01090.4 

07:57:06.64 

09:56:34.90 

2009-10-15 

2010-04-06 

1.33 

13.0 

B 



0805-077 

0808-0751 


08:08:15.54 

-07:51:09.89 

2009-10-22 

2010-04-12 

1.05 

12.7 

B 

... 


0804+499 

0808+4950 


08:08:39.67 

49:50:36.50 

2009-10-12 

2010-03-27 

1.30 

13.3 

B 

... 


0823+033 

0825+0309 

PKS0823+033 

08:25:50.30 

03:09:24.00 

2009-10-24 

2010-04-14 

1.48 

13.7 

B 

r- 


0827+243 

0830+2411 

OJ248 

08:30:52.08 

24:10:60.00 

2009-10-21 

2010-04-09 

1.48 


... 

... 


0836+710 

0841+7053 

4C71.07 

08:41:24.37 

70:53:42.20 

2009-10-12 

2010-03-17 

1.80 

12.9 

B 

20.5 

B 

0851+202 

0854+2006 

OJ287 

08:54:48.80 

20:06:30.00 

2009-10-27 

2010-04-15 

3.90 

13.1 

B 

21.8 

B 

0906+430 

0909+4253 

3C216 

09:09:33.50 

42:53:46.08 

2009-10-23 

2010-04-09 

1.02 

13.2 

B 

•«. 


0917+449 

0920+4441 


09:20:58.46 

44:41:53.99 

2009-10-25 

2010-04-11 

1.27 

12.7 

A 

21.7 

A 

0923+392 

0927+3902 

4C39.25 

09:27:03.01 

39:02:20.90 

2009-10-28 

2010-04-15 

7.06 


... 

22.7 

B 

0945+408 

0948+4039 

4C40.24 

09:48:55.34 

40:39:44.60 

zuuy-ju-3i 

2U1U-U4-1H 

1.32 

13.1 

B 

... 


0953+254 

0956+2515 


09:56:49.88 

25:15:16.05 

2009-11-08 

2010-04-27 

1.13 

13.1 

B 

... 


0954+658 

0958+6533 

S40954+65 

09:58:47.20 

65:33:54.00 

2009-10-21 

2010-03-30 

1.03 

12.8 

B 

21.7 

B 

* 1036+054 

1038+0512 


10:38:46.78 

05:12:29.09 

2009-11-29 

2010-05-16 

1.13 
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Table 1. continued. 


(1) 

Name 

(2) 

J2000 Name 

(3) 

Alias 

(4) 

R.A.(J2000) 

(5) 

Dec(J2000) 

(6) 

1st scan 

(7) 

2nd scan 

(8) 

[Jy] 

(9) 

logv. 

(10) 

Syn quality 

do 

logv,c 

(12) 

IC quality 

TEX 1040 +244 

1043+2408 


10:43:09.00 

24:08:35.00 

2009-11-19 

2010-05-10 

1.02 


ait 



1055+018 

1058+0133 

OL093 

10:58:29.61 

01:33:58.80 

2009-12-07 

2010-05-24 

3.96 

12.7 

B 

22.5 

A 

J1 130+3815 

1130+3815 

QSOB1 128+385? 

11:30:53.28 

38:15:18.50 

2009-11-21 

2010-05-15 

1.16 

12.9 

B 

... 


1150+812 

1153+8058 


11:53:12.50 

80:58:29.15 

2009-10-16 

2010-03-07 

1.11 



... 


1150+497 

1153+4931 


11:53:24.47 

49:31:08.83 

2009-11-16 

2010-05-11 

1.40 



... 

... 

1156+295 

1159+2914 

4C29.45 

11:59:31.83 

29:14:44.00 

2009-12-04 

2010-05-31 

2.00 

13.4 

A 

23.5 

B 

1219+044 

1222+0413 


12:22:22.55 

04:13:15.78 

2009-12-29 

... 

1.39 





1222+216 

1224+2122 

PKS 1222+2 16 

12:24:54.51 

21:22:47.00 

2009-12-18 


1.16 

14.2 

B 

22.7 

B 

1226+023 

1229+0203 

3C273 

12:29:06.69 

02:03:08.60 

2010-01-02 


24.37 

14.1 

B 

20.8 

A 

1228+126 

1230+1223 

3C274 

12:30:49.42 

12:23:28.00 

2009-12-26 


14.04 

13.4 

B 


»t« 

1253-055 

1256-0547 

3C279 

12:56:11.17 

-05:47:21.50 

2010-01-13 


17.20 

12.8 

B 

22.7 

B 

1308+326 

1310+3220 

AUCVn 

13:10:28.6 6 

32:20:43.80 

2009-12-21 


2.26 

12.6 

B 

... 


1324+224 

1327+2210 


13:27:00.86 

22:10:50.16 

2010-01-05 


1.15 

12.8 

B 


H , 

1413+135 

1415+1320 


14:15:58.80 

13:20:24.00 

2010-01-25 


1.51 


«*. 



1418+546 

1419+5423 

OQ530 

14:19:46.60 

54:23:14.00 

2009-12-06 


1.16 

... 

... 

20.6 

B 

1502+106 

1504+1029 

OR103 

15:04:24.98 

10:29:39.00 

2010-02-08 


1.30 





1510-089 

1512-0905 

PKS15I0-089 

15:12:50.53 

-09:05:59.00 

2010-02-14 

... 

2.48 

13*4 

B 



1546+027 

1549+0237 


15:49:29.44 

02:37:01.20 

2010-02-21 


2.42 

13.8 

B 

, 


1548+056 

1550+0527 


15:50:35.27 

05:27:10.45 

2010-02-21 


2.04 


... 



1606+106 

1608+1029 

4C 10.45 

16:08:46.20 

10:29:07.80 

2010-02-24 


1.65 

... 




1611+343 

1613+3412 

DA406 

16:13:41.00 

34:12:48.00 

2010-02-23 


3.03 

13.5 

B 

, 


1633+382 

1635+3808 

4C38.41 

16:35:15.49 

38:08:04.50 

2010-03-02 


2.90 

... 

... 

21.6 

B 

1638+398 

1640+3946 


16:40:29.63 

39:46:46.03 

2010-03-04 


1.18 


... 

«a» 


1642+690 

1642+6856 


16:42:07.85 

68:56:39.76 

2009-11-06 

2010-01-17 

1.74 



)T . 

, ,, 

1641+399 

1642+3948 

3C345 

16:42:58.81 

39:48:37.00 

2010-03-05 


7.89 

12.2 

A 

... 


1652+398 

1653+3945 

MARK501 

16:53:52.20 

39:45:36.00 

2010-03-09 


1.00 





1739+522 

1740+5211 

S41739+52 

17:40:36.98 

52:11:43.00 

2010-03-30 

... 

1.45 





1741-038 

1743-0350 

PKS 174 1-038 

17:43:58.86 

-03:50:04.60 

2009-09-13 

2010-03-18 

4.39 


T+ 

... 

... 

1749+096 

1751+0939 

PKS 1749+096 

17:51:32.70 

09:39:01.00 

2009-09-13 

2010-03-21 

4.50 

13.0 

A 

23.2 

A 

1803+784 

1800+7828 

S5 1803+784 

18:00:45.40 

78:28:04.00 

2009-10-11 

2010-02-02 

1.90 

13.5 

A 

23.1 

A 

1807+698 

1806+6949 

3C371.0 

18:06:50.70 

69:49:28.00 

2009-10-21 

2010-01-07 

1.37 

14.3 

A 

23.1 

B 

1823+568 

1824+5651 

4C56.27 

18:24:07.07 

56:51:01.50 

2010-04-17 


1.60 

13.0 

B 



1828+487 

1829+4844 

3C380 

18:29:31.80 

48:44:46.62 

2010-04-11 


2.38 

14.0 

B 



J 1849 15 +67064 

1849+6705 


18:49:15.89 

67:06:40.90 

2009-10-12 

2009-12-25 

1.82 

12.6 

A 


... 

1928+738 

1927+7358 

4C73.18 

19:27:48.50 

73:58:01.60 

2009-09-28 

2010-01-18 

2.61 



21.0 

B 

1954+513 

1955+5131 


19:55:42.74 

51:31:48.55 

2009-10-09 

2010-05-06 

1.35 

13.7 

B 

... 

... 

2007+776 

2005+7752 

S52007+77 

20:05:31.10 

77:52:43.00 

2009-09-25 

2010-01-27 

1.40 

13.1 

B 



2005+403 

2007+4029 


20:07:44.94 

40:29:48.60 

2009-10-17 

2010-05-02 

2.04 

... 

*** 



2021+614 

2022+6136 

OW637 

20:22:06.68 

61:36:58.80 

2009-12-20 

2010-05-26 

1.25 

... 

... 


... 

2037+511 

2038+5119 


20:38:37.04 

51:19:12.66 

2009-11-25 

2010-05-18 

2.46 

13.2 

B 



2121+053 

2123+0535 


21:23:44.52 

05:35:22.09 

2009-1 1-02 

2010-05-08 

1.74 

... 

... 


... 

2131-021 

2134-0153 

PKS2 13 1-021 

21:34:10.31 

-01:53:17.24 

2009-11-02 

2010-05-09 

1.46 





2134+004 

2136+0041 

0X057 

21:36:38.59 

00:41:54.21 

2009-11-04 

2010-05-10 

3.34 


t) 

21.0 

A 

2136+141 

2139+1423 


21:39:01.31 

14:23:36.00 

2009-11-11 

2010-05-16 

1.29 

12.5 

B 


... 

2145+067 

2148+0657 


21:48:05.46 

06:57:38.60 

2009-11-10 

2010-05-16 

7.37 





2200+420 

2202+4216 

BLLAC 

22:02:43.30 

42:16:39.00 

2009-12-13 

2010-06-06 

3.35 

118 

B 

22.6 

A 

2201+315 

2203+3145 

4C31.63 

22:03:14.98 

31:45:38.30 

2009-12-02 

2010-05-30 

2.53 

14.2 

B 

20.7 

B 

2201+171 

2203+1725 


22:03:26.89 

17:25:48.20 

2009-11-21 

2010-05-24 

1.21 

13.5 

B 
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Table 1. continued. 


0) 

Name 

(2) 

J2000 Name 

0) 

Alias 

(*) 

R.A.(J2000) 

(5) 

Dec(J2000) 

(6) 

1st scan 

(7) 

2nd scan 

(8) 

s™ [Jy] 

(9) 

log y, 

(10) 

Syn quality 

(ii) 

logvic 

(12) 

1C quality 

2216-038 

2218-0335 


22:18:52.04 

-03:35:36.90 

2009-11-14 

2010-05-21 

1.41 

14.0 

B 

... 

Mt 

2223-052 

2225-0457 

3C446 

22:25:47.26 

-04:57:01.40 

2009-11-15 

2010-05-23 

5.58 

12.7 

B 

22.0 

A 

2227-088 

2229-0832 


22:29:40.09 

-08:32:54.50 

2009-11-14 

2010-05-23 

1.94 

... 

... 

20.7 

B 

2230+114 

2232+1143 

CTA102 

22:32:36.41 

11:43:50.90 

2009-1 1-27 

2010-05-31 

3.84 

13.2 

B 

22.6 

B 

2234+282 

2236+2828 


22:36:22.47 

28:28:57.41 

2009-12-10 


1.03 

13.0 

B 

... 

,,, 

2251+158 

2253+1608 

3C454.3 

22:53:57.75 

16:08:53.60 

2009-12-06 

, ,, 

9.42 

13.9 

B 

22.6 

A 

4C45.51 

2354+4553 


23:54:21.68 

45:53:04.00 

2010-01-11 

... 

1.12 

... 


... 


2353+816 

2356+8152 


23:56:22.79 

81:52:52.26 

2009-09-19 

2010-02-08 

1.59 


... 
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Table 4. Simultaneity of high energy data to Planck observations displayed in the SEDs. 


B 1950 Name 

J2000 Name 

uvot 

SWIFT-XRT 

FERMI 

0003-066 

0006-0623 

sim.l st 

sim.lst 

sim_lst 

0007+106 

0010+1058 



27m 

0048-097 

0050-0929 

sim-1 st 

sim.lst 

sim.lst 

0059+581 

0102+5824 

m 

♦ »» 

sim-lst 

0106+013 

0108+0135 

»p* 


sim_lst 

JO 125— 0005 

0125-0005 



sim_lst 

0133+476 

0136+4751 

sim.lst 

sim.lst 

sim_lst 

0149+218 

0152+2207 


... 

sim-lst 

0202+149 

0204+1514 

~ri 


sirrulst 

0212+735 

0217+7349 

... 


27m 

0224+671 

0228+6721 

,,, 

... 

sim_lst 

0234+285 

0237+2848 

sim.lst 

sim.lst 

sim_lst 

0235+164 

0238+1636 

sim.lst 

sim.lst 

sim_lst 

0238-084 

0241-0815 

,,, 


sim_lst 

0306+102 

0309+1029 



27m 

0316+413 

0319+4130 



27m 

0333+321 

0336+3218 



27m 

0336-019 

0339-0146 



sim_lst 

0355+508 

0359+5057 



sim.lst 

0415+379 

0418+3801 



sim.lst 

0420-014 

0423-0120 

sim-lst 

sim.lst 

sim.lst 

0430+052 

0433+0521 

sim_lst 

sim.lst 

sim.l st 

0446+112 

0449+1121 

rn 


27m 

0458-020 

0501-0159 



27m 

0507+179 

0510+1800 



sim.l st 

0528+134 

0530+1331 

sim_2nd 

sim_lst_2nd 

sim_lst 

0552+398 

0555+3948 

sim_lst 

sim.lst 

sim.l st 

0605-085 

0607-0834 



sim-lst 

0642+449 

0646+4451 

sim_lst 

sim.lst 

sim.lst 

0716+714 

0721+7120 

sim.l st_2nd 

sim.lst.2nd 

sim.l st_2nd 

0723-008 

0725-0054 




0735+178 

0738+1742 



27m 

0736+017 

0739+0137 

sim_lst 

sim_lst 

sim-lst 

0748+126 

0750+1231 

sim_2nd 

sim_2nd 

2M_2nd 

0754+100 

0757+0956 



2M_lst 

0805-077 

0808-0751 



sim-lst 

0804+499 

0808+4950 

sim_lst 

sim_lst 

2M_lst 

0823+033 

0825+0309 

!>■ 


27m 

0827+243 

0830+2411 



2M_lst 

0836+710 

0841+7053 

sim_2nd 

sim_2nd 

2M.2nd 

0851+202 

0854+2006 

sim_2nd 

sim-2nd 

2M_lst_2nd 

0906+430 

0909+4253 




0917+449 

0920+4441 

sim.lst 

sim_lst 

2M_lst 

0923+392 

0927+3902 

sim_2nd 

sim_2nd 

sim_2nd 

0945+408 

0948+4039 



2M-1 st 

0953+254 

0956+2515 



2M_lst 

0954+658 

0958+6533 

sim.lst 

sim_lst 

2M_lst 

1036+054 

1038+0512 



2M.lst 

TEXl 040+244 1043+2408 



2M_lst 

1055+018 

1058+0133 

sim.lst 

sim.lst 

2M_lst 

J1 130+3815 

1130+3815 



2M_lst 

1150+812 

1153+8058 



sim-lst 

1150+497 

1153+4931 

sim.lst 

sim.lst 

2M_lst 

1 156+295 

1159+2914 

sim.l st_2nd 

sim_lst_2nd 

sim_lst_2nd 

1219+044 

1222+0413 



27m 

1222+216 

1224+2122 

sim.lst 

sim.lst 

sim-lst 

1226+023 

1229+0203 

sim.lst 

sim-lst 

sim-lst 

1228+126 

1230+1223 

sim.lst 

sim.lst 

sim_lst 

1253-055 

1256-0547 

sim.lst 

sim.lst 

sim-lst 

1308+326 

1310+3220 

sim.lst 

sim_lst 

sim_lst 

1324+224 

1327+2210 

4*+ 


sim_lst 

1413+135 

1415+1320 



sim.lst 

1418+546 

1419+5423 

TT 

sim-lst 

27m 

1502+106 

1504+1029 



27m 

1510-089 

1512-0905 

sim.lst 


sim.lst 

1546+027 

1549+0237 

sim_lst 

sim-lst 

sim-lst 

1548+056 

1550+0527 

... 


sim_lst 
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Table 4 . continued. 


B 1950 Name 

J2000 Name 

UVOT 

SWIFT-XRT 

FERMI 
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sim.lst 
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1613+3412 
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1635+3808 


sim_lst 

sim_lst 

1638+398 
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i , 
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0003-066 J0006-0623 0003-066 J0006-0623 




Fig. 18. Left-panel: The SED of the source 0003-066. Grey circles show the historical data. The red circles show data simultaneous 
to the Planck observations. The dotted and dashed lines show the second and third degree polynomials, respectively, fitted to the 
synchrotron and IC bumps in the SED. Right-panel: The radio spectrum of 0003-066. Red circles, LF data simultaneous to Planck . ; 
red stars, ERCSC LFI data; blue circles, HF data simultaneous to Planet, blue stars, ERCSC HFI data. The dashed and dotted lines 
are fits to simultaneous LF and HF data. 
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Fig. 20. 0048-097 
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Fig. 23. JO 125-0005 
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Fig. 24. 0133+476 
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Fig. 27. 0212+735 
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Fig. 33. 0316+413 
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Fig. 36. 0355+508 
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Fig. 44. 0552+398 
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Fig. 55. 0823+033 
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Fig. 56. 0827+243 
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Fig. 62. 0945+408 
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Fig. 68. J1 130+38 15 
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Fig. 69. 1150+812 
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Fig. 71. 1156+295 
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Fig. 74. 1226+023 
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Fig. 75. 1228+126 
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Fig. 77. 1308+326 
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Fig- 87. 1633+382 


1638+398 J 1640+3946 1638+398 J1640+3946 




Fig. 88. 1638+398 


1642+690 J 1642+6856 


1642+690 


J1 642+6856 



10.0 



10 l ° 10 " 
v (Hz) 


Fig. 89. 1642+690 


i l 1 1 ill 


Planck Collaboration SEDs and radio spectra of northern AGN, Online Material p 31 


1641+399 J 1642+3948 1641+399 J1642+3948 




Fig. 90. 1641+399 
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Fig. 96. 1807+698 

1823+568 J 1824+5651 
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Fig. 98. 1828+487 
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Fig. 99. J184915+670 
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Fig. 100 . 1928+738 
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Fig. 103 . 2005+403 
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Fig. 104 . 2021+614 
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Fig. 105 . 2037+511 
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Fig.107. 2131-021 
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Fig. 109 . 2136+141 
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Fig. 111. 2200+420 
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Fig. 112. 2201+315 
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Fig. 114. 2216-038 
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Fig. 115. 2223-052 
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Fig. 116. 2227-088 
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Fig.117. 2230+114 
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Fig. 118 . 2234+282 
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Fig. 120. 4C45.51 
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Fig. 121. 2353+816 


